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ABSTRACT 

Some aspects of gravitational lensing by large scale structure (LSS) are 
. investigated. We show that lensing causes the damping tail of the cosmic 
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microwave background (CMB) power spectrum to fall less rapidly with 
decreasing angular scale than previously expected. This is due to a transfer of 
power from larger to smaller angular scales which produces a fractional change 
in power spectrum that increases rapidly beyond t ~ 2000. We also find that 
lensing produces a nonzero mean magnification of structures on surfaces of 
constant redshift if weighted by area on the sky. This is a result of the fact 
that light-rays that are evenly distributed on the sky oversample overdense 
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regions. However this mean magnification has a negligible affect on the CMB 
power spectrum. A new expression for the lensed power spectrum is derived 
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and it is found that future precision observations the high-£ tail of the power 
'. spectrum will need to take into account lensing when determining cosmological 

parameters. 
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c/2 . 1. Introduction 

Previous discussions of gravitational lensing by large-scale structure have concentrated 
on calculating the shear and convergence along unperturbed light-paths, i.e. what the 
geodesies would be were there no fluctuations (e.g. |(Seljak 1992)| , and references cited 



therein). Three basic methods have been adopted. The first is by numerical simulation 
(e.g. |(Fukushige et al. 1994 )|) . This method often suffers from limited resolution and 



overly idealized cosmological models. Another method has been to use a model where light 
travels freely in a constant background density between clumps of localized mass densities 
(Fukushige et al. 1994 )| , |(Bessett et al. 1994)| . This is not considered to be a realistic 



cosmological model, because of the wide range of length scales on which galaxy clustering 
is observed. What appears to be the best method thus far is to take a smooth field of 
density fluctuations and calculate the shear and convergence along unperturbed light-paths. 
This can be done with the use of optical scalars |(Gunn 1997)| , |(Blandford et al. 1991)| or 



equivalently by using methods based on those of |(Kaiser 1992) 



In particular, |(Seljak 1992)| has applied the techniques of Kaiser |( Kaiser 1992 )| to 



the lensing of the CMB. He found that lensing results in a relatively small smoothing of 
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the CMB power spectrum which makes peaks and troughs somewhat less distinct. This 
smoothing is due to fluctuations in the magnification of structures on the surface of last 
scattering. The average magnification was assumed to be zero, as it is to first order. Seljak 
also found that evolving the deflecting density fluctuations by linear or nonlinear theory 
makes little difference in the results for I < 1000. 

We show here that deviations of the light-paths from their form in an unperturbed 
universe result not only in fluctuations in the magnification around a mean of zero, but also 
a shift in the mean to a positive value. Light-paths are attracted by regions of overdensity 
and repelled by regions of underdensity. This means that the column density of mass 
seen by the observer is larger on average than what would be expected using unperturbed 
light-paths. The predominantly positive second derivatives of the potential in overdense 
regions produces a shear between light paths which acts to magnify images. At the same 
time, the average shear between light-paths is, to a lesser extent, reduced by the increase 
in the density of light paths in overdense regions. The net result is that objects on surfaces 
of equal redshift or cosmological time will on average appear larger than in an unperturbed 
universe. The apparent violation of flux conservation can be resolved by realizing that the 
area of a surface of constant redshift is smaller when light-paths are perturbed. In angular 
size coordinates, light travels "slower" in regions of low potential. 

The other and more important aim of this paper is to show that after lensing the CMB 
power spectrum will be enhanced over the unlensed power spectrum at small angular scales 
or large I. Power is transferred upward in i in the damping tail. This result is independent 
of the existence of a nonzero mean magnification. The paper is organized as follows: In 
the next section we introduce the formalism used to calculate the lensing effects of LSS. In 
section 3 it is shown how lensing will change a generic CMB power spectrum. In section 4 
the formalism is applied to some specific cosmological models. 



2. Calculating the Magnification 

Throughout this paper, the Universe is assumed to have Robertson- Walker geometry 
together with small fluctuations. This implies that the density fluctuations are isotropic 
and the universe is homogeneous on the average. We also assume that the lensing is weak 
so that there are not multiple images of a single source. It can be shown without difficulty 
that the cross-section of regions with densities over the critical density required to produce 
multiple images is rather small, so that they should not play an important role in the 
statistical properties of lensing over large regions of sky (Schneider, Ehlcrs fc Falco 1992) 
(Kochanek 19951 . 



In the longitudinal gauge with conformal time, the metric takes the form 



ds 2 = a(r) 2 (-(l + 2(j))dT 2 + (1 - 2(P)dx 2 ); dx 2 = dr 2 + g{r) 2 (d 2 6 + sin 2 (#)d 2 0), (1) 
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where <p <C 1 and g(r) = { R sinh (r /R), r, R sin (r/R)} for the open, flat and closed global 
geometries respectively. The curvature scale is R = \H Q y/l — Q — ^a| _1 . Because Maxwell's 
equations are conformally invariant, for the purpose of finding light-paths the expansion of 
the universe can be ignored as long as conformal time is used. In general, light follows a 
geodesic that is a solution to: 

d dx v (\) 1 dx a (X) dxP(X) n dr . , x dt 

5A g "-ar = 2**- ix ix • p = Tx = ^ TX (2) 

Choosing A = r by normalizing p° and taking the unperturbed path to be the r-axis, the 
evolution equation to first order in the potential becomes 

<h<W -2g(r)- 2 <P,e = -2g(rr l d ± . (3) 



dr 2 

Since r = — r + r Q to first order in this equation can be solved as a function of r: 

S$ t ( r ) = 5 -^- = ^- f dr'g{r - r')<P A {r'). (4) 



g(r) g(r) Jo 

This must be evaluated along the path that the light bundle has followed. The first order 
effects arise from evaluating it along the unperturbed path. To find the correction due to 
the perturbation of the path, we expand the potential to first order: 



Sdi (r) = 5 -^- = f dr'g(r - /) U^r') + 5x k {r')<P, k {r ! )) (5) 

g(r) g[r) Jo v ' 

where the potential is now evaluated along the unperturbed path. Repeated indices are 
summed over the two components perpendicular to this path. Likewise the Sx(r') inside 
the integral can be approximated by the first order deflection calculated from equation (f|) 
evaluated along the unperturbed path. The shear tensor which measures the distortion and 
expansion of an infinitesimally thin beam is then 

= -qqJ = % + A % = ^r)l dr '9(r')9(r ~ r')<j> ^{r') (6) 

+JL f dr ' f dr " g {r _ r i) g (r> - r ") W)^")*^) + gir")^")^')] . 
g(r) Jo Jo 

In general this expansion is not justified for fluctuations of all scales. However it can be 
shown by explicit calculation that higher order terms are quite small in realistic models. If 
we assume that the relevant scales are much smaller than the curvature scale we can Fourier 
decompose the potential, 

0(r) = / t^0(A;, T = T - r ) e -*te e -i(rfcr+ ff (r)9.fcx). ( 7 ) 
J (27r) d 

In this section we assume that the angles involved are small enough that a local Cartesian 
coordinate system, 8, can be set up with the usual inner product. 
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The average value of the shear tensor can be found by substituting equation (^) into 
equation (|]) and using the assumption that the Fourier components are uncorrelated, i.e. 
{4>(k, r)4>(k', t)*) = (27r) 3 P^(/c, T)5(k — k'). In the case of the linear evolution of the potential 
fluctuations in a universe dominated by nonrelativistic matter, the time dependence of the 
potentials can be factored out of its Fourier components, <j>(k, r) = D(r)<fi(k). In this case, 

(<%) 

W(r, rk r , g(r)9 ■ k±) 



Equation (0) can be interpreted as consisting of two contributions. The term with g(r') 
is due to the average potential, as sampled by the light paths, being below average. The 
g(r") term results from the density of light paths being higher in areas of low potential. 
The reduced separation between light paths makes them converge less rapidly. The second 
term almost cancels the first term because in popular models the k values that contribute 
most are large enough that the oscillations of the exponential restrict r' — r" to be small. It 
appears that the coherence length of structure is small enough to make this magnification 
negligible. 

The time enters into these calculations because it is a function of the radial coordinate 
that parameterizes the light path. All the significant quantities calculated in this section, 
such as the second term in equation @, contain two integrations over this parameter. 
However, when \r' — r"|(= \t' — t"\) is large, larger then some ill-defined "coherence length", 
the potential fluctuations at these two points are uncorrelated and do not contribute 
significantly to the integrals. If the potential changes slowly enough it will not change 
significantly in the time it takes light to travel one "coherence length" and we can take 
{<j>(k,T)<f){k' ,t')*) = {2ir) 3 D(r = (t + t')/2) 2 P^(k)5(k- k'). We will call this the average 
time assumption. It can be avoided at the expense of complicating the evaluation of the 
integrals.^ These complications of time evolution are largely avoided in the flat-CDM model 
because the potential is time-independent in linear theory. 

The quantity of interest for applications to the CMB is the difference in the deflections 
of light paths that are observed to have an angular separation of s on our sky. This can be 
found by integrating the shear tensor 

= f /2 >\>,,iO)<io r (io) 

J— s/2 



1 If a significant amount of hot dark matter exists or there has been substantial nonlinear evolution, the 
factorization of <f> into k- and r-dependent parts will not be possible. In these models the whole power 
spectrum must be kept within the r integrals which in general must be evaluated numerically. 
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kikiP^Wir, rk r , g(r)6 ■ k A 



(2tt) 3 

dr' I dr"D(r')D(r" 



(9) 



xg(r - r')g(r' - r")\g(r') - gtf^e^^'^+bW-gV'W-kd^ 
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/3(s) has both a small average due to second order terms and a variance which is dominated 
by first order terms. Combining equations ([It]), (0) and (|6]), 

(P\\As) 2 ) = f (w)) 2 Jo dr' & dr" f ^klP^TMr - r')g(r - r »)e^r>-r») (n) 

x {J [w-} - J [w + ] ± (Ji[w_]/w_ - Ji[w + ]/w+)} , 

where to_ = k^sg(r' — r")/2 and u> + = k±sg(r' + r")/2. The plus sign in the second 
line is for /3m, the component parallel to s, and the minus is for the perpendicular 
component, (3±. In these coordinates the cross terms vanish. To reduce the numerical work 
necessary to integrate these oscillatory integrands it is useful to make an approximation. 
This approximation can be understood by first changing variables from {r', r"} to 
{r = (r' + r") /2, y = r' — r"}. When y is large the fluctuations in <fr will add incoherently. 
We assume that the coherence length is small enough that g(r — r') ~ g(r — r") ~ g(r — r) 
and that J {k±sy /2) ~ 1 for all relevant k± and s. Then the y integration of the exponential 
can be done giving 2rj Q (rk r ). This spherical Bessel function suppresses the contribution of 
modes with k r much greater than 1/r. It is usually the case in lensing situations that the 
peak in P^k) is at a k ^> 1/r so that everything but j (rk) can be brought out of the k r 
integral and k± replaced with k. In other words, the modes that contribute most are nearly 
perpendicular to the line of sight. The result is 

2 \ 2 r . f dk. Q 



P l± (s) 2 )^l^—j dr J —k*P^T)g{r -r) 2 {l - J [ksg(r)}± J x [ksg{r))/ksg{r)} . 

(12) 



This is equivalent to result that is found in |(Seljak 1992)| using the Fourier space 



Limber's equation derived by Kaiser |(Kaiser 1992)| . We find explicitly that for popular 



models, equation flT2"|) estimates equation (|TT| ) quite well for both large and small s, and 
appears to do well for intermediate values. 



3. The CMB Power Spectrum 

We wish to calculate the effect gravitational lensing should have on fluctuations in the 
CMB. Lensing will not change the temperature or surface brightness of the CMB, but it will 
change the size and shape of features. In this paper we restrict ourselves to the CMB power 
spectrum. In |(Seljak 1992)| a relation between the lensed and unlensed power spectra was 
derived by replacing the spherical harmonic transformation of the CMB fluctuations with 
its Fourier transform. This is a good approximation if I ^> 1 and Is C 1 for all the relevant 
scales of s and t. In the appendix we derive a relation that avoids the second assumption 
by remaining in spherical harmonic space. We also retain the anisotropic contributions that 
are dropped in |(Seljak 1992)| . We find that for the models tested the two methods agree 



very well if the anisotropic contributions are included in both cases. Given that the two 
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methods of calculation give the same results we prefer ours because we find it to be the 
faster one in our implementations. 

The transformation of the power spectrum derived in the appendix is: 

2f + 1 r 



Cf ~ yCe-—^\ rf S sin( S )P4cos( S )]{e- £ ' 2 ^^ 2 >/ 2 P^[x] (13) 
tt^q 2 Jo 1 



l'=0 

+- 



2 

where C^ b is the observed power spectrum. 



x=cos(s) 



4. Applications 

For the purpose of making some quantitative predictions, we adopt a flat, A = 0, cold 
dark matter cosmological model with adiabatic initial density perturbations. We use the 
CDM linear power spectrum given by |(Bardeen et at. 1986)| with the correction for finite 
baryon density given by |(Sugiyama 1995)|: 



P(k) = Ak n T(ke 2nb /h 2 ) 2 (14) 

2.34 9 



T(q) = ln(1 + 2 . 34t?) [1 + 3.89g + (16.1g) 2 + (5.46g) 3 + (6.71g)T 1/4 (15) 



This is related to the potential power spectrum by Poisson's equation. The 
Hubble parameter, H a = lOO/i km s _1 Mpc _1 and the COBE normalization is 
A = 6.6 x 10 5 /i~( 3+n ) Mpc 4 K Bennett et al. 1996)|. The nonlinear evolution of the 



power spectrum is calculated using the fitting formulae of |( Peacock fc Dodds 1996 



The second moments of (3\\(s) and /3±(s) are given in figure (|1|) along with an example of 
the first moment. It can be seen that at small angles (f3(s) 2 ) is approximately proportional 
to s 2 . At large separations the deflections of the two paths are uncorrelated, and so 
P\\(s) ~ P±(s) ~ (8(6) 2 ) which is independent of s. This constant could also be calculated 
using equation (U). In general increasing the Hubble parameter increases lensing effects 
because it shifts power down to smaller physical scales. 

Figure (|2]) shows the effect lensing has on the CMB power spectrum calculated using 
equation (0). In doing this calculation for large I one must take care that the Legendre 
polynomials and the integration are being accurately calculated. The unlensed spectra are 
produced using Seljak & Zaldarriaga's code, |(Seljak fc Zaldarriaga 199671 , up to ^ = 3000. 



Beyond this point a fitting formula for the damping tail t{l + l)Ci oc f^ 1 ) exp[-2(£/£ D ) m ] 
is used where m and £d are given in |(Hu fc White 1 996)1 The plots of Cf'/Ci shows not only 



that the fractional changes in the C/s are substantial at large £, but also that the damping 
tail of the observed and the intrinsic, unlensed spectra will have different functional forms. 
The calculation is also done with the J\ terms dropped from ([12]) so that the lensing is 



0.01 0.10 1.00 10.00 100.00 1000.00 

s(orcmin) 



Fig. 1. — The first two moments of f3(s) for some CDM models are given here. The straight 
line is ((3\\(s)) in a COBE normalized model. The two curves of each type are (j3\\(s)) (the larger) 
and (P±(s)) (the smaller). The COBE normalized models have only linear evolution of the matter 
power spectrum. All the models have h = 0.6 and ag is the rms density fluctuation in a sphere of 
radius 8/i _1 Mpc. 




1000 1500 2000 2500 3000 JBOO 1500 2000 2500 3000 3500 



Fig. 2. — The lensed CMB power spectrum for a CDM model with h = 0.6, fi = 1 and Q b = 0.04: 
The top panels are the lensed and unlensed spectra and on the bottom are their ratios. The left 
panels show n = 1.0, ag = 0.6 normalized models. The linear matter power spectrum evolution 
and the isotropic lensing approximations are shown. The right panels show the effect of tilting the 
power spectrum with COBE normalization and linear evolution. In the top right panel the solid 
curves are the unlensed spectra. 
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isotropic, (P\\(s) 2 ) = (j3±(s) 2 ). The results reproduces the ones we get from equation (A8) 
(Seljak 1992)| . When n, the spectral index of the primordial power spectrum, is reduced 



in 



the steepening of the damping tail partially compensates for the weaker lensing and results 
in a comparatively small change in C° b /C e ll. Including nonlinear structure formation is 
found to increase the lensing effect significantly for £ ^ 2000. The range of £-space that 
contributes in ([13]) is surprisingly large. At £ = 3000 about 30% of the change in power 
is due to contributions from £ < 1000 and 8% is from £ < 500, in the a 8 = 0.6 model. In 
this sense lensing can be thought of as transferring power from the acoustic peaks to the 
damping tail in such a way that the variance is conserved. 

The parameters m and £e> are dependent on h 2 Q , h 2 Qb, h 2 Q\, and the average 
temperature. Hu & White |(Hu fc White 199"6~)| have proposed that measuring m and 
£d would be a good way of determining these density parameters. However lensing will 
cause the spectrum to fall less rapidly than expected. If we fit a damping curve to the 
1000 < i < 3000 part of the as = 0.6-nonlinear lensed spectrum in figure (H) m and Id are 
changed by 19% and —17% from the unlensed spectrum. If measured in this way lensing 
will cause the angular size of the damping scale at the time of decoupling, ip , to be 
overestimated. Likewise a small value for m would mean that the thickness of the surface 
of last scattering would be overestimated. The change in m is especially significant since it 
is otherwise quite a weak function of cosmological parameters. Within the acceptable range 
of parameter space m changes by only about 10% in the unlensed spectrum. 



5. Conclusion 

We have shown that gravitational lensing can have a significant effect on the CMB 
power spectrum at small angular scales. The mean magnification will probably be too 
small to detect, but variations in the magnification will cause the damping tail to decrease 
less rapidly with increasing £. We have also found that nonlinear structure formation 
and anisotropic contributions the transformation of the power spectrum are important 
at large £. Acoustic peaks at large £ may be smoothed to such an extent that they are 
unidentifiable. The effects of lensing can be removed from the spectrum, but a model for 
both the lensing potential and the unlensed CMB power spectrum must be assumed. In 
addition the transformation of the power spectrum is nonlinear although it seems well 
behaved. This increases the amount of potential information in the damping tail, but makes 
the interpretation of future small-scale observations more ambiguous. Perhaps the power 
spectrum of density perturbations will be more tightly constrained by other means. 

In this paper we have only shown results for flat cosmological models with no 
cosmological constant, A. Lensing effects will be somewhat smaller in both low density 
and A models, because of the appearance of the mass density in Poisson's equation, 
P^ik.r) = 9a(r)~ 2 l^H„£;~ 4 P(A;, r)/4. This factor overcompensates for the increase in path 
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length and growth in fluctuations with lookback time (Scljak 1992) . With reasonable values 



for Q a and A the change in the spectrum is still significant. 

Interferometers are under construction that will be capable of probing the predicted 
CMB fluctuations from 150 to 3500 in t. The window size in £-space for an interferometer 
is ~ 2nD where D is the diameter of the dishes in units of the wavelength. The proposed 
instruments will operate at around 30GHz and have dish diameters of tens of centimeters so 
an £-space resolution of 100 should be achievable. This should allow for many independent 
measurements of the rate at which the tail falls with £. In addition, mosaicing over the 
sky can further narrow the window. These experiments will use multiple frequencies so 
that foregrounds such as the Sunyaev-Zel'dovich effect can be removed. We have shown 
that lensing corrections increase the amplitude by a factor of ~ 2 or more above £ ~ 3000 
in flat CDM models with Hubble constants in the observed range. With £-space windows 
in the above range any experiment that is capable of detecting the unlensed spectrum at 
these high £'s will be measurably affected by lensing. It will then be essential to include the 
lensing contribution to the CMB fluctuations in order to utilize the tail beyond £ ~ 2000 
for cosmological parameter estimation. 
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A. Appendix 

We start by expressing the observed temperature fluctuations, I o (0) = (T(8) — (T))/ (T), 
in terms of the spherical harmonic expansion of the intrinsic, unlensed fluctuations, 

oo I 

WA)=Y. E a lm YP{0'A') (Al) 

£=0m=-£ 

where {9', <fr'} is the position from which the light that is observed at position {9, <p} would 
have come from if there were no lensing. The two-point correlation function at an angular 
lag of s can be found with the help of (a^ m o|/ m /) = ^w^mm'^t an d the addition theorem for 
spherical harmonics, 

00 (2£ + 11 - 
C(s) = (I (9 l , <h)I o {0 2 , <h)) = E { — L2 C £ (P e {cos(\s + p(s)\)]), (A2) 

where the \s + (3(s)\ is the angle that would separate the two points on the surface of last 
scattering if there where no lensing. It has been assumed that the intrinsic fluctuations 
in the CMB are uncorrelated with the fluctuations in the potential that contribute to the 



-10- 



lensing. This is a very good approximation in realistic models. Because j3(0) = 0, the 
variance of the CMB is unchanged by lensing, i.e. the quantity J2eLo(2£ + 1)Q is conserved. 

Now the correlation function can be transformed back into the observed power spectrum 
using C° b = (27r) Jq ds sin(s)£(s)P^[cos(s)]. To express the C% b in terms of moments of (3(s) 
instead of the average of Legendre functions we expand these functions in a power series 

(P,[cos(|a + ?(*)!)]) * E ^^y^ ^ll^^M (A3) 



n=0 

1 

+ 2 



COS is 



mis?) -—(Msf) 



p'M- 



where the prime and (n) superscript represent the first and nth derivative with respect to 
x = cos(s). In each n-term, except n — 1, only the terms with the lowest powers of /3(s) 
are kept. Since the nth derivative is of order I larger than the (n — l)th derivative and I 
is large, the higher order terms are not important. Because the component of /3(s) parallel 
to s, (3\\(s) and the component perpendicular to s, (3±(s), are uncorrelated their cross terms 
are of order (ft(s) 2 ) 2n for small s and do not contribute significantly. For n = 1, higher 
order terms are kept because in practice they contribute significantly to the transformation 
of the power spectrum. 

Only the first term in the (3±(s) has been kept. A noticeable improvement in accuracy 
can be made by including higher order terms in the /3n(s) series. This is because at some 
point i becomes of order /3(s) _1 and Pe[x] oscillates with a period of order (3{s). This 
problem can be effectively circumvented if \ 0\\(s) — (/3||(s))J is a Gaussian random variable. 
Since /3||(s) is the result of a radial integral over many coherence lengths the central 
limit theorem supports the assumption of Gaussianity even if the potential field is weakly 
non-Gaussian. In practice {/3(s)) is small enough that we need only keep its first term in 
equation (|A5|). For the other terms Gaussianity requires (/3(s) 2m ) ~ (2m)!(/5(s) 2 ) m /2 m m! 
and for large £, P^ 2m \x] ~ (— l) m £ 2m P([x\/ {1 — x 2 ) m which can be easily shown directly 
from Legendre's equation. As will be seen the lensing has significant effects only at large I. 
With these simplifications 

Cf ~ rC ? — rf S sin( S )P4cos( S )]{e- £ ' 2 ^W 2 >/ 2 P,4x] (A4) 

77^n 2 JO l 



l'=Q 
1 

+ 2 



COS 



s)(m 2 )-—^?) 



p'AA-(P\\(s))sHs)p'i> 




It is clear from this expression that the change in the power spectrum will become significant 
when £ 2 (j3(s) 2 ) is significant. Neglecting the first moment of /3(s) and taking the small angle 
limit of the P^[x] term, which is a good approximation in practice, gives equation (|13|). 
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